First-principles study of ground state properties and high pressure behavior of Th02 
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The mechanical properties, electronic structure and phonon dispersion of ground state Th02 as 
well as the structure behavior up to 240 GPa are studied by using first-principles density-functional 
theory. Our calculated elastic constants indicate that both the ground state fluorite structure 
and high pressure cotunnite structure of Th02 are mechanically stable. The bulk modulus, shear 
modulus, and Young's modulus of cotunnite Th02 are aU smaller by approximately 25% compared 
with those of fluorite Th02. The Poisson's ratios of both structures are approximately equal to 0.3 
and the hardness of fluorite Th02 is 27.33 GPa. The electronic structure and bonding nature of 
fluorite Th02 are fully analyzed, which show that the Th-O bond displays a mixed ionic/covalent 
character. The valence of Th and O ions in fluorite Th02 can be represented as Th'^ ®^'*"^ and 
qO.452- rpj^g phase transition from the fluorite to cotunnite structure is calculated to be at the 
pressure of 26.5 GPa, consistent with recent experimental measurement by Idiri et al. For the 
cotunnite phase it is further predicted that an isostructural transition takes place in the pressure 
region of 80 to 130 GPa. 

PACS numbers: 61.50.Ah, 61.50.Ks, 71.15.Mb, 63.20.dk 
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I. INTRODUCTION 

Besides uranium and plutonium, thorium is also one 
kind of important nuclear materials. Metal thorium and 
its compounds have been widely investigated both experi- 
mentally and theoretically since metal thorium was found 
in 1828. Among thorium compounds, thorium dioxide 
Th02 , which is a stable diamagnetic transparent insulat- 
ing solid, has attracted special attention. In addition to 
its usage as an important nuclear fuel material, thorium 
dioxide has also been used as a solid-state electrolyte. 
In particular, due to its prominent hardness, Th02 has 
potential interests as an optical component material and 
laser host. 

Recently, there has occurred in the literature a series 
of experimental reports on pressure-induced phase tran- 
sition of Th02 P,[2,[l. At ambient pressure, Th02 crys- 
tallizes in the (CaF2-type) fluorite structure with space 
group Fm3m (No. 225). By using the energy dispersive 
x-ray diffraction (EDXRD) method, Dancausse et al. Q 
reported that at 40 GPa, Th02 undergo a phase transi- 
tion to an orthorhombic structure of cotunnite (PbCl2) 
type with space group Pnma (No. 62). Later, through 
improving experimental measurement technique, Idiri et 
al. [l| observed that this phase transition really begins 
around 30 GPa, and the two phases coexist in a wide pres- 
sure range of nearly 20 GPa. Additional high-pressure 
Raman spectroscopy measured by Jayaraman et al. [2\ 
also suggested that Th02 starts to transform into the 
cotunnite structure around 30 GPa. 

In contrast to the above-mentioned extensive experi- 
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mental investigation, a systematic theoretical study on 
the phase transition of Th02 at high pressure is still 
lacking. In particular, considering the obvious discrep- 
ancy between two experimental groups P, Q involving 
the transition pressure, such a theoretical investigation 
from basis quantum mechanics is not only complemen- 
tary but also indispensable. Motivated by this obser- 
vation, in this paper, we present a first-principles study 
by calculating the total energies and enthalpies of Th02 
at its experimentally established crystalline phases. Our 
calculation shows that the transition pressure is around 
30 GPa, consistent with the recent experiment by Idiri 
et al. [1]. 

The other task for this paper is to theoretically present 
a thorough description of physical, mechanical and chem- 
ical bonding properties of Th02 at its ground state of 
the fluorite phase. To date, theoretical studies of the 
ground-state behavior of Th02 are very scarce in litera- 
ture 0, 0, S 0, IB and some of them are even inconsistent 
with the experimental data to a large extent. For ex- 
ample, without taking into account the 5/ state, Kelly 
et al. [4] calculated the bulk modulus B of Th02 to be 
290 GPa, which is far from the experimental value of 198 
GPa. On the other hand, based on a purely ionic model, 
Harding et al. Q obtained 5=175 GPa, which, compared 
to the experimental data, clearly underestimate the bind- 
ing interaction in the material to a large extent. This as 
whole encourages us to theoretically report a systematic 
investigation on the ground-state properties of Th02 at 
its ambient phase. Besides the consistency with the ex- 
isting experimental data, we expect that the new results 
for the first time predicted in this paper and experimen- 
tally inaccessible at present, such as the hardness, the 
phonon spectrum, and the charge transfer of the Th-0 
bond, will be of great help for a primary understanding 
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of Th02. The rest of this paper is arranged as follows. 
In Sec. II the computational method is briefly described. 
In Sec. Ill and Sec. IV we present our calculated results, 
some of which are made comparison with available exper- 
imental and theoretical results. In Sec. V we summarize 
the conclusions of this work. 



II. COMPUTATIONAL METHOD 

Our total energy calculations are carried out by em- 
ploying the plane-wave basis pseudopotential method 
as implemented in Vienna ab initio simulation package 
(VASP) 0]. The exchange and correlation effects are 
described by the DFT within the generalized gradient 
approximation (GGA) [lOj. The electron and core in- 
teractions are included using the frozen-core projected 
augmented wave (PAW) approach which combines the 
accuracy of augmented-plane-wave methods with the ef- 
ficiency of the pseudo-potential approach [llj . The tho- 
rium 6s^7s^6p^6d^5f^ and the oxygen 2s^2p* electrons 
are treated as valence electrons. Note that although the 
5/ state are empty in elemental Th, this level turns to 
evolve into a hybridization with O 2p state in the va- 
lence band, as well as to prominently contributes to the 
conduction band (see Fig. 3 below). We use a 9x9x9 fc- 
point grid with Monkhorst-Pack scheme [l^ for fluorite 
Th02 and 9x15x9 fc-point grid for cotunnite structure. 
Electron wave function is expanded in plane waves up to 
a cutoff energy of 500 eV, and all atoms are fully relaxed 
until the quantum mechanical forces become less than 
0.02 eV/A. 

To obtain optimized lattice constants of ground state 
Th02, first, we calculate several total energies at dif- 
ferent lattice constants around the experimental value. 
Then we fit the energy-volume data to the third-order 
Birch-Murnaghan equation of state (EOS) to give the 
theoretical equilibrium volume, minimum energy, bulk 
modulus B, and pressure derivative of the bulk modulus 
B'. On the other hand, the bulk modulus B, shear mod- 
ulus G, Young's modulus E, Poisson's ratio v, can also 
be derived from the elastic constants. We find that the 
bulk modulus B obtained by these these two ways are 
in good agreement, indicating that our calculations are 
self-consistent. 



III. GROUND STATE PROPERTIES FOR 
FLUORITE-TYPE THORIUM DIOXIDE 

A. Atomic structure and mechanical properties 

At ambient condition thorium dioxide crystallize in a 
CaF2-like ionic structure. Its cubic unit cell is composed 
of four Th02 formula units with the thorium atoms and 
the oxygen atoms in 4a and 8c sites, respectively [see 
Fig. 1(a)]. Each Th atom is surrounded by eight O 
atoms forming a cube and each O connects with four 




FIG. 1: Cubic unit cell for Th02 in space group Fm3m (a) 
and orthorhombic unit cell in space group Pnma (b), larger 
black spheres stand for Th atoms and the smaller white O. 

Th atoms to build a tetragon. A particular feature of 
this structure is the presence of a large octahedral hole 
sited at position (i, i, i). The present optimized lattice 
constant (ao) is 5.62 A(see Table I), in good agreement 
with the experimental [HQ value of 5.6 A. 

Elastic constants can measure the resistance and me- 
chanical features of crystal to external stress or pressure, 
thus describing the stability of crystals against elastic 
deformation. For small strains Hooke's law is valid and 
the crystal energy i? is a quadratic function of strain 
[l3|.Thus, to obtain the total minimum energy for cal- 
culating the elastic constants to second order, a crystal 
is strained and all the internal parameters relaxed. Con- 
sider a symmetric 3x3 nonrotating strain tensor e which 
has matrix elements Sij {i,j — 1, 2, and 3) defined by 




Such a strain transforms the three lattice vectors defining 
the unstrained Bravais lattice {a^, fc = 1, 2, and 3} to the 
strained vectors ^15] {a'^.} as defined by 

al. = (I + £)a,., (2) 

where I is a unit 3x3 matrix. Each lattice vector afe or 
aj, is a 3 X 1 matrix. The change in total energy due to 
above strain (1) is [l3| 

l^E = E{V,{e,})~E{V,G) 

= vj^a^e., + I ^ a.e.e, H- ©({e^}), (3) 

2—1 — ^ 

where E{V, 0) is the total energy for the unstrained crys- 
tal, (Ti are the elements of stress tensor, and are the 
elastic constants. For cubic structure of Th02, there are 
three independent elastic constants, i.e., Cn, C12, and 
C44, which are calculated through a proper choice of the 
set of strains {ei,i = 1,...,6} listed in Table II. In our 
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TABLE L Calculated lattice parameter, elastic constants, bulk modulus B, pressure derivative of the bulk modulus B , shear 
modulus G, Young's modulus E, Poisson's ratio v, hardness, and energy band gap (Eg) for cubic Th02. As a comparison, other 
theoretical works and available experimental data are listed. 



Property 


This study 


Previous calculation 


Exp er iment 




VASP-GGA 


FP-LTMO-GGA 




nn (A) 


5.62 


5.61" 




Cii (GPa) 


349.5 


376.0" 


367' 


Ci2 (GPa) 


111.4 


109.8" 


106"^ 


C44 (GPa) 


70.6 


68.1" 


79^= 


B (GPa) 


191 


198", 198^* 


198(2)\ 195(2)'' 


b' 


4.5 


4.5", 4.2"^ 


4.6(3)\ 5.4(2)'* 


G (GPa) 


87.1 


94.1" 


97'= 


E (GPa) 


226.8 


243.8" 


249'= 


V 


0.302 


0.295" 


0.285" 


Hardness (GPa) 


27.33 






Eg (eV) 


4.673 


4.522"^ 


6^ 



" Reference |8|], Reference [1, " Reference pj], '' Reference @|, '^Reference 



TABLE IL Three strain combinations in the strain tensor 
[Eq. (1)] to calculate the three elastic constants of cubic 
ThOa. 



Strain 


Parameters (unlisted e,;=0) 


AE/V in 0(S'^) 


£3 


ei=5 62 —5 63=5 

6i=5 63=5 
64=5 65=5 66=5 


f(Cii +2Ci2)5" 
(Cn + Ci2)52 

|C445^ 



first-principles calculations, the strain amplitude S is var- 
ied in steps of 0.006 from (5=-0.036 to 0.036 and the total 
energies E{V, 6) at these strain steps are calculated, and 
then fitted through the strains with the corresponding 
parabolic equations of AE/V as given in Tables II to 
yield the required second-order elastic constants. Note 
that while computing these energies all atoms are allowed 
to relax with the cell shape and volume fixed by the 
choice of strains {cj}. Our calculated elastic constants 
for the ground-state fluorite structure of Th02 are listed 
in Table HI As a comparison, previous theoretical re- 
sults based on the full-potential linear mufhn-tin orbital 
(FP-LMTO) method [S], as well as available experimen- 
tal data[l^ are also presented in Table I. It is clear that 
our calculated elastic constants are in good agreement 
with the experimental results. 

After obtaining elastic constants, we can calculate bulk 
and shear moduli from the Voigt-Reuss-Hill (VRH) ap- 
proximations [17, 18', 19]. The Voigt bounds [13, Ho] on 
the effective bulk modulus By and shear modulus Gy 
are 



By = (Cii+2Ci2)/3 



(4) 



and 



Gv = (Cn - C12 + 3C44)/5. (5) 
Under Reuss approximation [18], the Reuss bulk modulus 



Bji and Reuss shear mudulus Gn are 
Bb = By 

and 



(6) 



Gr = 5(Cn - Gi2)Gm/[4Cm + 3(Cii - C12)]. (7) 

The bulk modulus B and shear modulus G, based on 
Hill approximation (Toj , are arithmetic average of Voigt 
and Reuss elastic modulus, i.e., B=^{Bji + By) and 
G—^{Gr + Gy). The Young's modulus E and Poisson's 
ratio V for an isotropic material are given bv[2l| 



and 



WG 

3B + G' 



3B - 2G 
2{3B + G)' 



(8) 



(9) 



The calculated results for these moduli and Poisson's 
ratio for the fluorite Th02 are listed in Table I. Note that 
we have also calculated the bulk modulus B by fitting the 
Murnaghan equation of state. The derived bulk modulus 
turns out to be exactly the same as that from the above 
VRH approximation, which again indicates that our cal- 
culations are consistent and reliable. For comparison, 
other theoretical results d, Q and available experimental 
datad @, [l6j are also shown in Table I. It can be seen 
that on the whole, our present results compare well with 
those previous experimental and FP-LMTO theoretical 
results. Concerning the Poisson's ratio, it is well known 
that for the common materials that have much smaller 
shear moduli compared with the bulk moduli, their Pois- 
son's ratio is close to 1/3. For the present Th02 system, 
our calculated shear modulus G is much lower than the 
bulk modulus G. Thus, our calculated result of 0.302 
for the Poisson's ratio, as well as the derived result of 
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0.285 according to the experimentally determined elastic 
constants 16] and using Eq. (9), can be well understand- 
able. 

Furthermore, hardness is also one fundamental physi- 
cal quantity when considering the phase stability and me- 
chanical properties. Note that the hardness is important 
for the applications of Th02 as both nuclear material and 
optical components. So here we calculate the hardness 
of Th02 by using the approach raised by Simunek et al. 
[23 |. In the case of two atoms 1 and 2 forming one bond 
of strength si2 in a unit cell of volume f2, the expression 
for hardness has the form [l^l 

H = {C/n)bi2Si2e^''f\ (10) 

where 



Si2 = V (eie2)/(nin2di2), Ci = Zi/ri (11) 

and 

/2 = C-^^f = 1 - [VT^/(ei + 62)]^ (12) 

ei + 62 

are the strength and ionicity of the chemical bond, re- 
spectively, and di2 is the interatomic distance; C=1550 
and CT=4 are constants. The radius is chosen to make 
sure that the sphere centered at atoms i in a crystal 
comtains exactly the valence electronic charge Zi. For 
fluorite structure Th02, fei2=32 counts the interatomic 
bonds between atoms Th (1) and O (2) in the unit cell, 
711=8 and n2=4 are coordination numbers of atom Th 
and O, respectively, ri = 1.81 (A) and r'2=1.00 (A) are the 
atomic radii for Th and O atoms, respectively, Zi — 12 
and Z2—6 are valence charge for Th and O atoms, re- 
spectively, c?i2=2.43 (A) is the interatomic distance, and 
57=145.53 (A^) is the volume of unit cell. Using Eqs. 
(10)-(12), we obtain Si2=0.081 and /2=0.0025. The hard- 
ness of Th02 at its ground-state fluorite structure is thus 
given by H=27.3 (GPa). This indicates that the fluorite 
Th02 is a hard material and approaches to a superhard 
material (hardness > 40 GPa). The high hardness of this 
crystal can be understood from the dense crystal struc- 
ture, which results in high valence electron density and 
short bond distances. The unusual mixture of covalent 
and ionic components in the Th-0 chemical bond will be 
discussed in the following subsection. 



B. Electronic structure and charge distribution 

Almost all the macroscopical properties of materials, 
such as hardness, elasticity, and conductivity, originate 
from their electronic structure properties as well as chem- 
ical bonding nature. Therefore, it is necessary to perform 
the electronic structure analysis of Th02. The calculated 
band structure (left panel) and total density of states 
(DOS, right panel) of fluorite Th02 are shown in Fig. 
[21 The present calculated energy band gap of 4.673 eV 
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FIG. 2: Band structure (left panel) and total density of states 
(DOS, right panel) for ground state Th02 with Fermi energy 
level Ef taken at eV as shown by the dashed lines. 



(listed in Table [T|, which is consistent with previous FT- 
LMTO theoretical result of 4.522 eV [Q], indicates that 
Th02 is a typical insulator. The underestimation of band 
gap compared with experimental value of 6 eV (23j is due 
to the drawback of the exchange-correlation approxima- 
tion (GGA). As a comparison, the orbital- resolved partial 
densities of states (PDOS) for one Th atom and two O 
atoms in a unit cell are also presented in Fig. [3l One can 
see that the low bands covering from -18.33 to -11.22 eV 
mainly consist of O 2s and Th 6p states which shows a 
clear hybridization. This hybridization is also responsible 
for the covalency of Th02 (see further discussion below). 
The high valence bands located just below the Fermi level 
are mainly contributed by O 2p states with a little Th 
6c? and 5/ states and have a band width of 3.79 eV. This 
illustrates a significant charge transfer from Th 6d and 
5/ states to O 2p states. As for the conduction bands, 
the DOS is mainly featured by Th 5/ states, mixed with 
a little Th 6d and O 2p states and has a width of 1.78 
eV. 

In order to gain more insight into the bonding nature of 
ground state Th02, we also investigate the valence charge 
density distribution. The calculated valence charge den- 
sity map of the (110) plane is plotted in Fig. |4l It is 
obvious from Fig. 4 that the charge density around Th 
and O ions are all near spherical distribution with slightly 
deformed toward the direction to their nearest neighbor- 
ing atoms. Furthermore, the charge density around Th 
and O ions is high while there is almost no valence charge 
in the large octahedral-hole interstitial region. This sug- 
gests that remarkable ionization exists for thorium and 
oxygen ions and significant insulating property exhibits 
in Th02. In order to describe the ionicity quantitatively 
and more clearly, we plot in the inset in Fig. [4] the line 
charge density distribution along the nearest Th-0 bond. 
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FIG. 3: (Color online) Partial density of states (PDOS) for 
ground state TI1O2 at GPa. The Fermi energy level is zero. 




FIG. 4: Valence charge density of Th02 in (110) plane. There 
▲ stands for Th, A for O atoms. The contour lines are drawn 
from 0.0 to 1.5 at 0.1 e/A'' intervals. The line charge density 
distribution between Th atom and the nearest neighbor O 
atom is shown in the inset, where the arrow indicates the 
minimum value. 



Therein the arrow indicates the charge density minimum 
of 0.45 e/A^. The distance between this charge den- 
sity minimum point and Th ion is 1.32 A. According to 
this, we can define the ionic radius of Th in Th02 as 
rTii=1.32 A. After subtracting this value from the Th-0 
bond length (2.43 A), one obtains the oxygen ionic ra- 
dius in Th02 as 7'o=l.ll A. The corresponding valency 
of Th and O in Th02 can be estimated by calculating the 
valence charges within the spheres of ionic radii. After 
integration, we find that there are 8.166 electrons around 
Th and 6.452 electrons around O. As a result, the valency 
of Th and O can be formally represented as Th^'*^"*"^ and 
qo.452-^ respectively. This indicates that the Th and O 
atoms in Th02 are ionized inequivalently. In other words, 
Th is ionized almost as Th^+, while O is weakly ionized. 




FIG. 5: Calculated phonon dispersion curves along the high- 
symmetry directions (left panel) and total density of states 
(DOS, right panel) for ground state Th02. 

This fact, on the other hand, suggests that besides ion- 
icity, the Th-0 bond also includes a prominent covalent 
components, which is consistent with the above PDOS 
analysis. In fact, the minimum value 0.45 e/A^ of charge 
density along the Th-0 bond in Th02 is prominently 
higher than that in typical ionic crystals. Based on the 
PDOS and charge density analysis, therefore, we derive 
that the Th-0 bond is a mixture of covalent and ionic 
components. From this aspect, it is now understandable 
that the previous full-ionic model Q has largely under- 
estimated the bulk modulus of Th02 . 



C. Phonon dispersion curve 

The calculation of phonon frequencies of the crystalline 
structure is one of the basic aspects when considering the 
phase stability, phase transformations, and thermody- 
namics of crystalline materials. By using the Hellmann- 
Feynman theorem and the direct method, we have calcu- 
lated the phonon curves along some high-symmetry direc- 
tions in the Brillouin zone, together with the phonon den- 
sity of states. For the phonon dispersion calculation, we 
use the 2x2x2 fee supercell containing 96 atoms and the 
4x4x4 Monkhorst-Pack fc-point mesh for the Brillouin 
zone integration. In order to calculate the Hellmann- 
Feynman forces, we displace four atoms (two Th and two 
O atoms) from their equilibrium positions and the am- 
plitude of all the displacements is 0.03 A. The calculated 
phonon dispersion curves along the T-X-W-L-T direc- 
tions are displayed in Fig. O The T-X, X-W, and L-T 
lines are along <0 1>, <1 2>, and <1 1 1> directions, 
respectively. 

For Th02 in the CaF2-structure primitive cell, there 
are only three atoms. Therefore, nine phonon modes ex- 
ist in the dispersion relations. One can see that there 
is no gap between the optic modes and the acoustic 
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branches and the LO-TO also has no sphtting at T point. 
Due to the fact that thorium is heavier than oxygen atom, 
the vibration frequency of thorium atom is lower than 
that of oxygen atom. Therefore, the phonon density of 
states can be viewed as two parts. One is the part lower 
than 5.5 THz where the main contribution comes from 
the thorium sublattice, while the other part higher than 
5.5 THz, included both acoustic and optical branches, are 
dominated by the dynamics of the light oxygen atoms. 

IV. HIGH PRESSURE BEHAVIOR OF THO2 

In order to investigate the high pressure behavior of 
Th02, we have optimized the structural parameters of 
its Fm3m and Pnma phases at different pressures by 
using GGA method. To aviod the Pulay stress prob- 
lem, we perform the structure relaxation calculations at 
fixed volumes rather than constant pressures. For FmSm 
phase, due to its high symmetry, the structure relaxation 
calculations are performed at fixed volumes with no re- 
laxation of coordinates. However, for Pnma phase, the 
coordinates of atoms and the cell shape are necessary to 
be optimized due to their internal degrees of freedom. 

The optimized structural lattice parameters a, b, and 
c for the Pnma phase at GPa are 6.174, 3.776, and 
7.161 A, respectively, giving 1^=167.0 A'^. This volume is 
prominently smaller than the equilibrium volume of 177.5 
A'^ for the Fm3m phase. Using the same method as that 
has been emp loyed in our previous work on orthorhom- 
bic BeH2 [2J| , we have applied the strains and calculated 
the elastic constants, various moduli, and Poisson's ra- 
tio V for cotunnite-type Th02 at GPa. The results 
are collected in Table IIIII From Table IIIII the following 
prominent features can be seen: (i) The orthorhombic 
Th02 at zero pressure is mechanically stable because its 
elastic constants satisfy the following mechanical stabil- 
ity criteria [25j . 

Cii > 0, C22 > 0, C33 > 0, C44 > 0, C55 > 0, Cgg > 0, 
[Cn + C22 + C33 + 2(Ci2 + Ci3 + C23)] > 0, 
(Cii + C22 - 2C12) > 0, (Cn + C33 - 2C13) > 0, 
(C22 + C33 - 2C23) > 0. (13) 

In fact, one can see from Table IIIII that the calculated 
C12, C23, and Ci3 are largely smaller than Cn, C22, 
and C33. Therefore, the mechanical stability criteria is 
easily satisfied; (ii) Although the equilibrium volume of 
cotunnite-type Th02 is distinctly compressed, its bulk 
modulus, shear modulus, and Young's modulus are all 
smaller by approximately 25% than those in fluorite-type 
Th02. This intriguing behavior has also been observed 
in the UO2 system 26], which, like Th02, will undergo 
a Fm,3m—^Pnma phase transition. Since the 5/ or- 
bital band is occupied in UO2, while it is almost empty 
in Th02, thus this similarity of the softening in mod- 
uli upon phase transition for the two systems is clearly 
unrelated with the 5/ orbitals. From this aspect, we 



speculate that most actinide dioxides share this phase- 
transition involved similarity. Again, we have confirmed 
for the Pnma phase that the bulk modulus B calculated 
by fitting the Murnaghan equation of state equals to that 
through VRH approximation; (iii) The Poisson's ratio v 
of Th02 in Pnma phase is nearly the same as in Fmim 
phase, i.e., close to 1/3. This is understandable since the 
shear modulus is much smaller than the bulk modulus in 
the both two phases. 

Now let us see the phase transition energetics of Th02. 
The total energies (per unit cell) of the two phases at dif- 
ferent volumes are calculated and shown in Fig. [6l Obvi- 
ously, the Fm3>m phase is stable under ambient pressure 
while under high pressure the Pnma phase becomes sta- 
ble. According to the rule of common tangent of two en- 
ergy curves, a phase transition at 26.5 GPa is predicted 
by the slope shown in the inset of Fig. [6l Besides, we also 
determine the phase transition pressure by comparing the 
Gibbs free energy as a function of pressure. At K, the 
Gibbs free energy is equal to the enthalpy H, expressed as 
H=E+PV. Figure [7] shows the relative enthalpies of the 
cotunnite phase phase with respect to the fluorite phase 
as a function of pressure. The crossing between the two 
enthalpy curves in Fig. [7] readily gives phase transition 
pressure of 26.5 GPa, which is fully consistent with the 
above result in terms of the common tangent rule. This 
value is well close to the experiment measurement by Idiri 
et al. [1] and and by Jayaraman et al. [2| as ~30 GPa. 
Here, the minor theoretical underestimation by ^--^5.0 GPa 
is speculated to be caused by the possible existence of an 
energy barrier with an amplitude of Aw. To overcome 
this energy barrier, the external pressure P' should be 
larger than the conventional common-tangent-rule deter- 
mined pressure P by an amount AP=P' —P. According 
to the experimentally measured transition pressure p' m 
and our theoretically obtained common tangent curve in 
Fig. [SI we deduce the energy barrier amplitude to be 
~0.06 eV (per formula unit). This value is too small to 
survive the Pnma phase to ambient conditions, which is 
consistent with the fact that to date, the Pnma phase of 
Th02 has only be observed under high pressures. 

Figure [5] shows the relative volume V/Vq evolution 
with pressure for Th02 in both Fm3m and Pnma phases. 
For comparison, the experimental data are also shown in 
the figure. Clearly, our calculated P-V equation of state 
is well consistent with the experimental measurement for 
the both two phases of Th02. Specially, at the calcu- 
lated transition pressure (26.5 GPa), our result in Fig. [8] 
gives that the volume collapse upon phase transition is 
5.9%. This value is very close to the experimental data of 
6.1% [11 with a little bit underestimation, which is obvi- 
ously within the accuracy of GGA. Note that unlike the 
cubic Fm3m phase, in which all coordinated are com- 
pletely determined at each volume due to the symmetry, 
the orthorhombic Pnma phase at each volume has two 
additional internal degrees of freedom that must be fully 
relaxed to obtain the energy minimum. Thus, we have 
further calculated and plotted in Fig. [5] the pressure de- 
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TABLE III: Calculated elastic constants, bulk modulus B, pressure derivative of the bulk modulus B , shear modulus G, Young's 
modulus E, Poisson's ratio v for cotunnite-type Th02 at GPa. Except the Poisson's ratio, all other values are in units of 
GPa. 
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FIG. 6: (Color online) Comparison of total energy vs the 
cell volume for Th02 in Fm3m and Pnma phases. A phase 
transition at 26.5 GPa is predicted by the slope of the common 
tangent rule, as shown in the inset. 
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FIG. 7: (Color online) Calculated enthalpy differences of 
Pnma phase with respect to FmZm phase as a function of 
pressure. 



pendence of the three lattice parameters (with respect to 
their equilibrium values) for the Pnma phase of Th02. 
In the experimentally reported pressure range from '^30 
GPa to '^80 GPa, our calculated evolution of the relative 
lattice parameters in the Pnma phase are well consistent 
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30 60 90 120 150 180 210 240 
Pressure (GPa) 

FIG. 8: (Color online) Calculated compression curves of Th02 
compared with experimental measurements. The volume col- 
lapses at our predicted phase transition point 26.5 GPa and 
experimental phase transition pressure 37 GPa are labeled. 



with the experimental observation [1]. In this pressure 
region, one can see from Fig. [3] that the responses of the 
three relative lattice parameters to the compression are 
anisotropic in the sense that the compression of the mid- 
dle axis a is most rapid compared to those of the long 
axis c and small axis &, which vary upon compression al- 
most in the same tendency. When the pressure becomes 
higher to be between 80 and 130 GPa, remarkably, it 
reveals in Fig. [9] that all the three relative lattice param- 
eters undergo dramatic variations by the fact that the 
small axis b has a strong rebound and the middle a is col- 
lapsed. When the pressure is beyond 130 GPa, then the 
variations of the three relative lattice parameters become 
smooth and approach isotropic compression. This signi- 
fies a typical isostructural transition for the Pnma phase 
of Th02. It should be stressed that this isostructural 
transition is not unique for Th02. Similar phenomenon 
has also been observed in other actinidc dioxides [2^ . 

To see the change of electronic structure of Th02 under 
high pressure, we present in Fig. [TU] the orbital-resolved 
PDOS for the Fm3m and Pnma phases at the theoretical 
transition pressure of 26.5 GPa. One can see that there 
occurs an understandable narrowing in the band gap 
(4.6 eV for the Fm3m phase and 3.15 eV for the Pnma 
phase). The valence and conduction bands carry a size- 
able downward shift to stabilize the high-pressure Pnma 
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Therefore, the transition-accompanied volume collapse in 
Th02 is irrelevant to the 5/ orbital effect. Second, In 
other actinide dioxides, such as Pu02, although taking 
into account the 5/ electronic interaction is crucial to 
obtain more reasonable ground state [2^, this local cor- 
relation effect is not expected to play a role in explaining 
the transition-accompanied volume collapse in these 5/- 
occupied actinide dioxides. The reason is simply because 
that inclusion of 5/ electronic localization will weaken 
the binding of cations and anions, causing an increase 
instead of a decrease in the volume of system. Therefore, 
the phenomenon of volume collapse during high-pressure 
phase transition of the actinide dioxides is mainly at- 
tributed to the ionic (instead of electronic) response to 
the external compression. 



FIG. 9: (Color online) Pressure behavior of the relative lattice 
parameters of the Pnma phase, where the drastic change in 
the relative lattice constants (region between dashed lines) 
indicates an isostructural transition. 
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FIG. 10: (Color online) Partial density of states (PDOS) for 
Fm3m phase (a) and Pnma phase (b) both at 26.5 GPa. The 
Fermi energy level is zero. The energy gaps of Fm3m phase 
and Pnma phase are 4.584 eV and 3.153 eV, respectively. 

phase. In addition, the valence and conduction bands 
are more widened and smooth in the Pnma phase. How- 
ever, these changes in the electronic structure are trivial 
and thus cannot be associated with the exotic transition- 
accompanied volume collapse in Th02 and other actinide 
dioxides. Concerning the 5/ electron effect, here we 
would like to present our viewpoint as follows. First, the 
5/ orbital is almost empty for both two phases of Th02. 



V. CONCLUSION 

In summary, the ground state properties as well as the 
high pressure behavior of Th02 were studied by means 
of the first-principles DFT-GGA method. The elastic 
constants and their derived moduli and Poisson's ratio 
were calculated for both the ambient Fm3m and the high- 
pressure Pnma phases of Th02 and were shown to accord 
well with recent reliable experiments. Based on these re- 
sults, mechanical stability of the Pnma phase at zero 
pressure was predicted. The hardness, lattice vibration 
dynamics, and the Th-0 chemical bond of the ambient 
phase were calculated and analyzed in order to support 
the practical application of of Th02. We showed that 
the Th-0 bond displays a mixed ionic/covalent character, 
with the valence of Th and O ions represented as Th'^-®^'*"'' 
and 0'^ *^^^. Here the ionicity is mainly featured by 
charge transfer from Th 6d/5f states to O 2p states, while 
the covalency is manifested by hybridization of oxygen 2 s 
and thorium 6p states. This mixed ionic/covalent feature 
makes Th02 a hard material, with its hardness calculated 
in this paper to be ~27 GPa. As another main task, we 
studied phase transition of Th02 under high pressure. 
Our calculated FmSm^Pnma transition pressure is 26.5 
GPa, according well with recent experimental results of 
~30 GPa. Our calculated transition-accompanied vol- 
ume collapse of 5.9% is also in good agreement with the 
experimental data of 6.1%. Under more higher pressure, 
we further found that there occurred an isostructural 
transition between 80 and 130 GPa for the Pnma phase, 
which is to be experimentally verified in future. 
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